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The pairing correlations in superfluid ''He are strongly modifled by quasiparticle scattering off 
a surface or an interface. We present theoretical results and predictions for the order parameter, 
the quasiparticle excitation spectrum and the free energy for thin films of superfluid ''He. Both 
specular and diffuse scattering by a substrate are considered, while the free surface is assumed to 
be a perfectly reflecting specular boundary. The results are based on self-consistent calculations of 
the order parameter and quasiparticle excitation spectrum at zero pressure. We obtain ne-w results 
for the phase diagram, free energy, entropy and specific heat of thin films of superfluid '^He. 



I. INTRODUCTION 

The superfluid phases of bulk '^He are condensates of 
p-"wave, spin-triplet Cooper pairs. The order parame- 
ters describing the A- and B-phases spontaneously break 
orbital- and spin-rotational symmetries, as ■well as dis- 
crete symmetries of space- and time-inversion, of the 
normal Fermi-liquid phase of '^Hei Pairing transitions 
of this type are called "unconventional" , and such sys- 
tems exhibit novel phenomena associated with the spon- 
taneously broken symmetries of the order parameter. A 
generic feature of superfluids with an unconventional or- 
der parameter is their sensitivity to scattering by impu- 
rities, defects and boundaries. Scattering of quasipar- 
ticles by these objects leads to the suppression of the 
order parameter (pairbreaking) , suppression of the su- 
perfluid transition as well as more subtle physical effects 
associated with the interplay of scattering and particle- 
hole coherence in the superfluid state. These effects 
persist over several coherence lengths, which for '^He is 
Co = hvf/2iTkBTc ~ 10^ - 10^ A. When ^He is placed in 
a container or confined to a geometry with dimensions of 
order the coherence length scale the effects of the bound- 
aries on the superfluid extend to all parts of the liquid, 
and there is in this sense no "bulk" phase. 

Superfluidity of '^He films was first reported in 1985>2, 
Basic properties of superfluid films, such as the tran- 
sition temperature, critical current, superfluid density, 
have been measured in several laboratoriesi^i^ Evidence 
of a presumed A-B phase transition has been reported 
by several groups, and a phase diagram has been con- 
structed over a limited range of temperatures and film 
thickness&i^ Experimental observation of third sound 
in ■^He films also shows anomalies in the mode spectrum 
as a function of temperature and film thickness which 
cannot be accounted for within the hydrodynamic theory 
applicable to superfiuid "^He^ There are many open ques- 
tions about the nature of the superfluid phases of '^He in 
restricted geometries, including the identiflcation of the 
superfluid order parameter and whether or not additional 
phases may be stabilized depending on the geometry and 
surface structure of the confining geometry. 

Theoretical models and calculations of the surface struc- 
ture and excitation spectrum of superfluid '^He films are 



important for understanding both thermodynamic and 
nonequilibrium properties of superfluid films. Detailed 
comparison between theory and experiment for the tem- 
perature dependence of the heat capacity and entropy, 
allows us to determine the density of states and low ly- 
ing excitations of the film. The heat capacity and en- 
tropy also enter the hydrodynamic equations that de- 
scribe damping of low frequency collective modes of su- 
perfluid films. 

In this paper we report a theoretical study of superfluid 
■^He in films with thickness ranging from D ^ \ — 15^0- 
This is a system similar to He in a slab, except the '^He is 
in general confined between different interfaces. Essential 
to any theory of superfluid ^He in conflned geometry are 
the boundary conditions that describe the effects of sur- 
face scattering on the pairing correlations, the quasipar- 
ticle spectrum and quasiparticle distribution functions in 
the case of nonequilibrium properties. 

Early theoretical investigations focused on pairbreaking 
of the order parameter near a wall in the Ginzburg- 
Landau (GL) regime, and the implications of the pair- 
breaking effect on the boundary conditions for the hy- 
drodynamic variables describing the A-phase, in partic- 
ular the £ vectoriiS These calculations, as well as calcu- 
lations of the suppression of the transition temperature 
for superfluid '^He in confined geometry^ were based 
on deGennes' formulation of inhomogeneous superfiuid- 
ity in terms of semiclassical correlation functions and a 
heuristic model of surface roughness which interpolated 
between specular and diffuse scattering of quasiparticles. 
More recent analyses based on the GL theory for the 
phase diagram and dynamical properties of '^He in slabs 
and cylindrical pores are described in Refs. ll^lT^ 

Extensions of surface pairbreaking calculations beyond 
the GL limit require a more detailed theoretical formula- 
tion of inhomogeneous states of superfiuid '^He. The most 
powerful theory of superfluid '^He is based on the qua- 
siclassical transport equationsii^ii^ti^ This theory is the 
natural extension of Landau's theory of normal Fermi liq- 
uids to include BCS pairing correlations. The quasiclas- 
sical theory is applicable to a broad range of phenomena 
and nonequilibrium states of inhomogeneous "^He. The 
central objects of the quasiclassical theory are the prop- 
agators that describe both the quasiparticle excitations 
of the condensed phases and the correlated pairs that 
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form the condensate. Theoretical calculations based on 
the quasiclassical theory for the surface order parameter 
and excitation spectrum near a wall or interface require 
boundary conditions for the quasiclassical propagators, 
generally formulated from scattering theory and a spe- 
cific model for the surface or interface. 
Boundary conditions describing reflection from an atom- 
ically rough surface were developed by Buchholtz and 
Raincr~ and implemented in the form of the randomly 
rippled wall (RRW) approximation for ^He-Bii2ii£ Alter- 
native formulations of diffuse boundary conditions were 
implemented by Zhang et ali^ based on scattering from a 
thin layer of atomic-size impurities coating an otherwise 
smooth surface, and by Thuneberg, et aliSi based on scat- 
tering from a distribution of "randomly oriented mirror" 
(ROM) surfaces. Boundary conditions describing rough 
surfaces which are neither perfectly specular nor fully 
diffuse were developed by Nagato et alm^ in terms of a 
random S-matrix. In the diffuse scattering limit the re- 
sults of Ref. [2^ for the order parameter suppression in 
^He-B agree well with those obtained by Zhang et al^ 
based on the TDL model. 

We report a theoretical analysis of the structure, exci- 
tation spectrum and thermodynamic properties of supe- 
fluid •^He films. The free surface is modelled as a spec- 
ular surface, and the film resides on a substrate that we 
assume is atomically rough. We model the scattering 
of quasiparticles by the substrate by introducing a thin 
layer of atomic-size impurities randomly distributed on 
the surface (TDL model). The width, d, of the impu- 
rity layer is assumed to be much less than the superfluid 
coherence length, while the mean free path Zjmp for quasi- 
particles propagating in the layer is much smaller than 
the width d. Thus, quasiparticles are strongly scattered 
inside the layer, but eventually scatter out of the layer at 
an angle uncorrelated with the incident trajectory. The 
limit d//imp ^ cxD as d — *■ describes a rough surface 
in the diffuse scattering limit. The specfic formulation 
of the impurity model for diffuse scattering that we use 
was introduced by Ovchinniko\i24 for diffuse scattering 
from atomically rough metallic interfaces, and was im- 
plemented for superfluid '^He by Kopnin^ 

We start from the quasiclassical theory of superfluidity 
with the boundary conditions described above and calcu- 
late the equilibrium properties of films of superfluid "^Hc 
in the weak-coupling limit. We calculate suppression of 
the order parameter in the film, the quasiparticle density 
of states, superfluid free energy, entropy and heat ca- 
pacity. We are especially interested in films of thickness 
D < I?AB, where £)ab ~ 9^0 marks the transition from 
the B-likc phase for thicker films to the A-like phase for 
thinner films. It is in this region that most film experi- 
ments have been done. Wc find that because of diffuse 
scattering at one of the interfaces, a band of subgap states 
is formed and these states exist throughout the film for 
sufficiently thin films. The thermodynamic properties of 
the A-phase are changed significantly by these gapless 
excitations. 



The main sections of this paper are organized as fol- 
lows. In section ^ we describe the theoretical formu- 
lation of the quasiclassical transport equations and the 
boundary conditions that we use in our calculations. Re- 
sults for the thin film phase diagram are summarized in 
section Unl In section Hvl we discuss the excitation spec- 
trum for both specular and diffuse scattering, while in 
section |V] we present the results for the free energy, en- 
tropy and heat capacity of thin films. More technical 
aspects of computing the free energy and implementing 
the diffuse boundary condition in the Riccati formulation 
of the transport equations are included in appendices. 

II. THEORETICAL MODEL 

We assume that the ^He film is on a substrate that is 
atomically rough on a scale much shorter than coherence 
length, ^0- The film has a well-defined thickness, D, that 
is larger than the atomic scale. The free surface of the 
film is assumed to be atomically smooth, c.f. Fig. ^ We 
assume there is negligible evaporation and vapor above 
the film. Thus, the liquid in the film is essentially under 
zero pressure. In this model '^He quasiparticles are spec- 
ularly reflected at the free surface. We also assume that 
film is invariant under translations and rotations in plane 
of the film (xy-plane) ^ and thus the physical properties 
of the film depend only on z-coordinate, which is normal 
to the substrate and the free surface. 




z = D 



z = 



FIG. 1: A thin film of ^He on an atomically rough substrate 
(at z = 0) with a free surface {z = D) which is specular. 

The calculations we report were carried out using the 
quasiclassical theory for superfluid •^He?^ supplemented 
by boundary conditions for surface scattering at the 
vapor-liquid interface and the film substrate. The cen- 
tral object of the quasiclassical theory is the propagator, 
g(p,R;£m), which is a 4 x 4 Nambu matrix - denoted 
by a widchat - in the combined particle-hole and spin 
spaces, and is defined in terms of an integration of the 
full Nambu propagator, 

G(p,R;e™) = j'^ dr e'''-'' J d^r e-'P '' x (1) 
- (T,«'(R + r/2,T)*(R-r/2,0)), 
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over a shell, \vf{p — pf)\ < Sc ^ Ef, in momentum space 
near the Fermi surface, 



?(p,R-; 



d£,v Tl3G(p, R; Sr. 



(2) 



The propagator is normalized by dividing by the weight 
of the quasiparticle pole in the spectral function, a. We 
use the Matsubara representation to calculate equilib- 
rium properties; the fermion Matsubara frequencies are 
Em = (2to + l)-KkBT. The four-component Nambu 
field operators are defined in terms of the bare fermion 
field operators hy — (-01, V'ij "^ji a-nd 'l'(r, r) = 
'I'1'(r,-T). 

For pure spin-triplet pairing the quasiclassical propa- 
gator g may be parametrized in particle-hole space by 
2x2 spin matrices for the diagonal (quasiparticle) and 
off-diagonal (Cooper pair) propagators, 



9 = 



g f 
f 9 



g + &-g {icrcFy) ■ f 
[iaya) • f 5 + o-"' • g 



(3) 



The 2x2 spin matrices are denoted by ordinary hats, e.g. 
g. The spin vector, a = (^yi ^z), are the Pauli matri- 
ces. One deviation from the matrix notation is p, which 
denotes a unit vector in the direction of the Fermi veloc- 
ity, v/(p) = Vfp. The components of the quasiclassical 
propagator are not all independent. The upper and lower 
particle-hole components are related by symmetries that 
follow from the fermion anticommutation relations. 



/(p, R; e™) = /(-p, R; e^)* = -/(p, R; -e™)^ , 
g{p, R; £,„) = g{-p, R; e,„)* = g(-p, R; -£,„)"' , 



(4) 



where 5'' is the matrix transpose of g. 
The quasiclassical transport equation that governs the 
evolution of the quasiclassical propagator, ^(p, R;£m), 
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[i£mT3 - A(p, R) , g(p, R; £,„)] 

+ iv/(p)- V5(p,R;£rn) =0, (5) 

with a constraint given by Eilenberger's normalization 
condition on ^, 



7(p,R;£,„)2 = -TT^l, 



(6) 



We have omitted the Landau molecular field self-energy 
in the transport equation, and we consider the pairing 
self-energy. A, in the weak-coupling limit, which is a con- 
venient choice for the order parameter. It is off-diagonal 
in particle-hole space, 



A 



A 
A 





iay& ■ A* 



iacTy ■ A 







(7) 



and parameterized by a spin-triplet order parameter de- 
fined by the vector, A(p,R). In the weak-coupling limit 



the order parameter is determined by the off-diagonal 
pair amplitude, f(p,R;£„i), from the gap equation, 

A(p, R) = T ' ^ / ^ V{p, p') f (p', R; £„) , (8) 

m 

where V{p,p') is the interaction in the spin-triplet pair- 
ing channel. For pure p-wave pairing we retain only the 
attractive £ = 1 interaction, V — 3Vip ■ p' . The cut-off, 
£c, and interaction, Vi, arc not measurable, but they are 
related to the bulk transition temperature by, 
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In 



1.13£c 



(9) 



which is used to eliminate the cut-off and pairing inter- 
action in favor of the measured bulk transition tempera- 
ture, Tc- 

The order parameter must be determined self- 
consistently with the solution of the transport equation 
for the propagator. This procedure, and the quasiclassi- 
cal transport equation, can be simplified by introducing 
a parametrization for the propagator that satisfies the 
normalization condition by construction and reduces the 
number of independent components. 
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1 + aa 



2a 



—2d —1 — aa 
where the prefactor is given by 
(1 - aa)-i 



N = 









(1 — aa)' 



(10) 



(11) 



The amplitudes d and a are 2x2 matrices in spin space 
which obey matrix Ricatti equations^i'SS'^S 



ivf ■ Vd + 2i£,„d — dAd -f A = 
zv f ■ Vd — 2i£„jd — dAd + A = . 



(12) 



We refer to d and d as the Ricatti amplitudes. The two 
Ricatti amplitudes are related to the particle- and hole- 
like projections of the off-diagonal propagators. 



-{in- g) V' , a= (^tt + g) V, 



(13) 



where the projection operators for the particle- (P+) and 
hole-like (P_ ) sectors are given by 



2'^ 
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For the case of spin-triplet pairing in zero field these am- 
plitudes can be parameterized as. 



d = (id' ■ Uy) ■ a and d = {icTy&) ■ a . 



(15) 



The Ricatti amplitudes are also related to each other by 
a symmetry that follows from symmetry relations for the 
propagators in Eqs. 10}, 



d(p, R; £„)* = d(-p, R; £„) . 



(16) 
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The Ricatti equations are easily integrated numerically, 
are numerically stable and provide a more efficient ap- 
proach to solving the quasiclassical transport equations 
than the "explosion method" Equations ()12(l are solved 
by integration along classical trajectories - forward for 
a and backward for a - starting from an initial value. 
The Ricatti equations must be supplemented by bound- 
ary conditions at the two interfaces. We do not have a 
bulk region in the "^Hc film, so generally we start from 
an arbitrary initial value at the free surface and compute 
along a classical trajectory with multiple reflections un- 
til the Ricatti amplitude at the surface converges. The 
integration procedure is described in more detail in App. 

HI 

The boundary conditions for the Ricatti amplitudes at 
the two interfaces are obtained from boundary conditions 
for the quasiclassical propagators. Specular reflection at 
the free surface requires matching of the propagators at 
the free surface for two trajectories, p and p, which arc 
related by p = p — 2fi(fi • p), 



5(p, £»;£,„) = dip, D; Sr. 



(17) 



Then the Ricatti amplitudes arc also matched at the sur- 
face in the same way. 



a(p, = a(p,D;e„ 



(18) 
(19) 



The boundary condition for the quasiclassical propaga- 
tor at an atomically rough surface is more complicated. 
A physical model for an atomically rough surface is pro- 
vided by a "thin dirty layer" (TDL) model for surface 
roughness obtained by coating a specular surface with 
a layer (of thickness d) of randomly distributed impuri- 
ties characterized by a mean free path of /jmpi^ In the 
TDL model the ratio, p = d/limp, describes the degree 
of surface roughness. For p — we recover a specularly 
reflecting surface, while p ^ oo corresponds to the fully 
diffuse surface. In the fully diffuse limit we implement 
Ovchinikov's boundary condition, which is a special case 
of the diffuse limit of the TDL boundary condition. The 
Ovchinnikov boundary condition requires self-consistent 
determination of the Green's function at the diffuse sur- 
face, g(p,0;£m). For outgoing trajectories (p^ > 0) the 
boundary condition for the Ricatti amplitude is. 



a(P,0) = -(iTT - ^tdl) Vt 



(20) 



where 5TDL(£m) and /TDL(£m) SuTC the propagators deep 
in the dirty layer, and which are related to the surface 
propagator by. 



E?TDL (^m) 



|p,|g(p,0;e„). (21) 



Pz >0 
Pz <0 



We give a short derivation of this boundary condition in 
appendix ^ 



Order Parameter 

We consider two possible phases in supcrfluid •^Hc films 
that have the same, or nearly the same symmetry, as 
the A- and B- phases of bulk supefluid '^He when one 
restricts the orbital symmetry group to SO{2). The order 
parameter for the B-like phase is of the form, 

Ab = (A|| (z)p„ A,| iz)py, A^{z)p,) , (22) 

while that for the A-likc ('axial') phase is given by, 



Aa = (0,0, A\\{z){p^+ipy)). 



(23) 



where || and _L refer to orbital motion, characterized by 
the direction of the relative momentum p, parallel and 
perpendicular to the surfaces of the film. The planar 
phase is a special case of the B-phase with A_l = 0. 
For the A-like phase we have A a || £ = z in order to 
minimize the nuclear dipolar energy. For the B-phase 
the order parameter in Eq. H22() is multiplied by a spin- 
orbit rotation matrix, TZ{n, t?), that is fixed by the dipole 
energy. 

The spatial profiles of the order parameter components 
are shown in Fig. [21 for both the B- (left panel) and A- 
like (right panel) phases. The dashed lines correspond 
to a film with two specular surfaces, while the solid lines 
represent a film with diffuse scattering from a substrate 
at z = and specular reflection from the free surface. 




diffuse and speculai' 
both surfaces specular 

J \ , I , 

4 6 



FIG. 2; The order parameters for the B-like and the A-like 
phases of ^He in a thin film. The coherence length, = 
hvf /2TikBTc, is approximately equal 73 nm at zero pressure. 

The orbital components of the order parameter that are 
perpendicular to the film interface, A_l, are suppressed at 
both interfaces. This suppression is related to the change 
of sign of Aj^ps when a quasiparticle is reflected by the 
surface. The parallel component, A||, is suppressed by 
diffuse scattering at the substrate in both phases. In the 
S-phase A^ is slightly increased for diffuse scattering 
because some of the spectral weight that is lost from Ay 
is transferred to A^. 
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The orbital structure of the order parameter for the A- 
phase leads to a simplification for the boundary condi- 
tion at the diffuse substrate. If we parameterize the off- 
diagonal component of g by 



(0, 0,/|i(p,)(p, + ip,))^ 



(24) 



then the angular integration in Eq. H21|l gives /tdl = 0. 
Thus, we have an explicit value for the Ricatti amplitude, 
d, at that surface, 

a(p,0;e™) = 0, p. > . (25) 

The quasiclassical Green's function at z = is then 



g(p,0,e„) 



1 

-2a -i 



P.>0. (26) 



The boundary value in Eq. H25|l speeds up numeri- 
cal integration since the calculation of a{pz > 0) and 
a(pz < 0) are now initial value problems; we start at 
z = with a{pz > 0, 0; Em) = 0, or a{pz < 0, 0; £,„) = 
and integrate the Ricatti equations directly to obtain 
a{pz > 0, z; Em) and a{pz < 0, z; Em)- 



III. 



PHASE DIAGRAM 



At zero pressure bulk '^He is in the superfluid B-phase for 
temperatures below Tc = 0.93 mK. When we confine the 
superfluid to a slab between two surfaces, or form a film 
on a substrate, we observe changes in the superfluid as we 
decrease the film thickness, D. The phase diagram of "^He 
in superfluid films, as far as it is known, is shown in Fig. 
13 Several phase transition lines calculated theoretically 
are shown, as well as points indicating possible phase 
transitions based on anomalies in several experiments. 

If one starts from the bulk superfluid B-phase and 
then reduces the film thickness, D, at constant tem- 
perature we expect to cross at least two phase bound- 
aries. As D is reduced the perpendicular component, 
A^(z), is suppressed, and at a critical film thickness, 
DAuiT), Aj_{z) vanishes. This signifies a transition to 
the planar phase with an order parameter of the form 
Ap = A\\{z){'Px ,Py , 0). The component A_l(z) van- 
ishes continuously, so this transition is second order. In 
the weak-coupling limit the planar and axial phases are 
degenerate. However, in bulk "^He strong-coupling cor- 
rections lower the free energy of the axial phase relative 
to the planar phase. Thus, if strong-coupling effects also 
stabilize the axial phase relative to the planar phase in 
a thin film then the second-order transition from the B- 
phase to the planar phase is pre-empted by a transition 
from the B-phase to the axial A-phase. Measurements of 
the heat capacity jump in bulk ■^Hc indicate that strong- 
coupling corrections are small at zero pressure, thus, the 
AB transition in the film is likely to be very weakly first 
order. With the exception of possible fine structure of 
the phase diagram close to the second-order transition 



line, Daj^{T), and properties such as the latent heat of 
transition, calculations of the thermodynamic properties 
of thin films based on the weak-coupling approximation 
are expected to be accurate. 

The perpendicular component of the order parameter is 
suppressed to zero even for a specular wall, so the AB 
transition occurs even in a film bounded by two specular 
surfaces. For a film on a rough substrate, the suppression 
of A|| by diffuse scattering leads to a small enhancement 
of A_l(z). As a result the B to A transition requires 
slightly thinner films for a rough substrate. This result, 
although the detailed shape of the phase boundary is 
slightly different, agrees with the calculations reported by 
Nagato et alm^ based on a different theoretical model for 
the surface roughness. However, NMR measurements^ on 
thin slabs of superfluid ^He show that the AB transition 
occurs at larger values of film thickness than predicted 
by the weak-coupling theory. This may indicate that the 
first-order AB phase boundary needs to be calculated 
with leading order strong-coupling corrections included 
in the theory, even at zero pressure. 




FIG. 3: Phase diagram for superfluid ^He films. The thick 
solid line represents the AB transition for a film in contact 
with a rough substrate and a specular free surface. The thin 
solid line is the AB phase boundary for a film with two spec- 
ular surfaces. The inset shows an enlarged portion of the 
AB phase boundary, where the second-order transition is re- 
entrant, A— >B— >A, as a function of temperature. The dashed 
line shows the normal (N) to A-phase boundary, with suppres- 
sion of the superfluid transition, T)!''" < Tc, resulting from 
diffuse scattering by the substrate. The individual points cor- 
respond to observed anomalies in measurements that may in- 
dicate a phase transition in the superfluid film: a) thin dashed 
line - anomaly in Ps{T, D)j~ b) open squares - mode anomaly 
in third soundf^ c) open circle - flow anomaly;- d) open dia- 
monds - flow anomalyi- 

An interesting feature of the calculated weak-coupling 
AB phase boundary below T/Tc w 0.4 is shown in the 
inset of Fig. O For films in contact with either a specular 
or rough substrate the second-order phase boundary is re- 
entrant as a function of temperature for a narrow range 
of film thicknesses. For example, for a film on a rough 
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surface with D ~ 9.4^0i upon decreasing the temperature 
below r^"'" the A to B transition occurs at Tab ~ O.SSTc. 
As the temperature drops further a re-entrant B to A 
transition occurs at Tba — 0.23Tc. Whether or not this 
re-entrance will survive strong-coupling corrections is not 
known. The re-entrance may also signal that a transla- 
tionally invariant A or B phase is unstable to formation 
of an inhomogeneous phase with lower free energy. In 
any event the fine structure of the phase diagram at low 
temperatures near D ~ 9.5^0j and the possibility of new 
phases stabilized by strong-coupling corrections, or which 
spontaneously break translation symmetry in the plane 
of the film, is outside the scope of this paper. 
For films that are thinner than 9.5 — 10 the planar 
or axial A-phase is the stable phase relative to the B- 
like phase. For the purpose of calculating the thermody- 
namic properties we assume that strong-coupling correc- 
tions stabilize the A-phase relative to planar phase in the 
film, however, this is really an open question. Strong- 
coupling corrections to the free energy for phases with 
strong spatial variations, as occurs in thin films, have not 
been calculated, so the relative stability of the planar and 
A-phase in thin films is unknown, either theoretically or 
experimentally. 

If the substrate were an ideal specularly reflecting sur- 
face then the superfluid A-phase would persist for film 
thicknesses approaching a few monolayers, or until the 
Fermi liquid properties and pairing interaction were mod- 
ified by finite size effects. But, the pair-breaking ef- 
fect of scattering off a rough substrate suppresses tran- 
sition temperature into the A-phase and at a film width 
Dan{T), which is substantially smaller than Dj^b{T) the 
superfluid A-phase is destroyed. The calculated transi- 
tion temperature for diffuse scattering is shown in Fig. 
13 This phase boundary was calculated by identifying the 
temperature and film thickness where the order parame- 
ter vanishes. We also obtained the superfluid transition 
from the calculated free energy by a least square fit of 
the known Ginzburg-Landau form for the free energy, 

a(r-r«'")2. 

Calculations of the transition temperature in thin slabs 
of ^He were carried out by Kjaldman et ali^ using a lin- 
earized gap equation and deGennes' formulation of the 
kernel in terms of the classical limit for the normal-state 
current-current correlation functioUfiE Our calculations 
agree well with the results for a slab if we take into ac- 
count that the width of the thin film of "^He-A is equiva- 
lent to a •^He-A slab of twice the width of the film. 

It should be noted that the phases considered here, even 
for thin ■^Hc films, assume thicknesses, D ^ A. We 
do not consider the two-dimensional limit of one or two 
atomic layers of ^He atoms on the surface of a substrate. 
The properties of 2D superfiuid ^He, if it exists, are ex- 
pected to be influenced by the reduced dimensionality. 
Ising-like, as well as Kosterlitz-Thouless-type transitions 
are predicted for 2D superfluid ^He-ASi^^. 

One additional note: observing the equilibrium phase 
boundaries may be complicated by metastability. Even 



though the planar phase, i.e. the B-phase with Aj^ = 0, 
and the axial A-phase are degenerate in weak-coupling 
they are unrelated by symmetry, and therefore, sepa- 
rated by an energy barrier. Thus, once established, the 
axial A-phase will be metastable with respect to the B- 
phase. The calculation of the barrier and the correspond- 
ing metastability lines in the phase diagram would pro- 
vide an important result, but are outside the scope of 
this article. 



IV. DENSITY OF STATES 

Pair-breaking by surface scattering leads to quasiparticle 
states below the gap. These excitations play an impor- 
tant role in the thermodynamic and transport properties 
of thin films of superfluid "^He. The sub-gap excitations 
are surface Andreev bound states. The mechanism lead- 
ing to their formation is closely related to the formation 
of bound states in the core of a vortexi2^ Andreev bound 
state formation occurs when the order parameter changes 
sign or phase along a quasiparticle trajectory. In the case 
of surface scattering the incident and reflected trajecto- 
ries generally correspond to very different order parame- 
ters; this is typically the case for an unconventional order 
parameter which breaks rotational symmetry. 

For example, consider the B-phase in contact with a 
specular surface. For an incident trajectory normal to the 
interface, p || z the B-phase order parameter changes sign 
upon reflection, i.e. A(p) = — A(p) for p ^ P = —p. 
This sign change leads to multiple Andreev reflection that 
generates a surface bound state at the Fermi level, i.e. 
s — 0. For trajectories away from normal incidence the 
components of the B-phase order parameter correspond- 
ing to orbital motion in the plane are present and do not 
change sign upon reflection. As a result the surface An- 
dreev bound state disperses as a function of the incident 
and reflected angles relative to the interface normal. 



FIG. 4: Skew scattering by the substrate. 

The bound state energy dispersion can be calculated 
approximately by neglecting the suppression of the or- 
der parameter at the surface, and assuming that surface 
scattering occurs on a cone defined by the angle 9 from 
the a;?;-plane, i.e. (0, ^ (tt — 6*, (/)2), = 02 ^ '/'i as 
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shown in Fig. 0] We then find bound-state poles in the 
retarded propagator, for either the B- or A-phase, given 

by 

£h = ±A|| sin 6* cos ^ . (27) 

The density of states (DOS) can be calculated once the 
order parameter and Landau molecular fields have been 
determined self-consistently. The most detailed informa- 
tion is contained in the angle-resolved local density of 
states, which is obtained from the diagonal component 
of the retarded quasiclassical propagator, 

N(p,R;e) = --amg^(p,R;e) (28) 

TT 

where (?^(p,R;e) is found by solving the quasiclassical 
transport equation for real energies; i.e. ie^ e -f- iO+ 
and the known order parameter and molecular fields. The 
local density of states for the B-phase near a wall shows 
quasiparticle states which develop below the bulk gap, 
and are bound to the surface, i.e. their spectral weight 
vanishes a few coherence lengths away from the surface. 
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FIG. 5; The angle-resolved local DOS for the ^He-B near a 
specular surface. The spectrum is calculated for T = O.STc. 
For clarity we have broadened the Andreev bound states with 
a width parameter, 77 = 10~''Ao. 

For example, the angle-resolved spectrum of supcrfiuid 
■^He-B near a specular surface, calculated numerically for 
a sclf-consistently determined order parameter, is shown 
in Fig. [S] For the specular reflection the position of 
the positive energy surface bound state depends on the 
angle of the incident trajectory, 9, approximately as Cb = 
A|| sinO. For normal incidence the bound state is at zero 
energy and disperses towards, eventually merging into, 
the continuum edge as the incident trajectory approaches 
grazing incidence. There is also a weak dispersion in 
the continuum edge reflecting the enhancement of Aj^ 
by surface scattering. 

At an atomically rough surface diffuse scattering couples 
an incident trajectory to all outgoing trajectories. This 



leads to mixing of states with different energies and thus 
to a band of sub-gap states for a given incident trajec- 
tory as shown in Fig. El The suppression of Ay for diffuse 
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FIG. 6: The angle-resolved local DOS for the ^He-B near 
a diffuse surface. The spectrum is calculated for T = O.STc. 
For clarity we have broadened the Andreev bound states near 
grazing incidence with a width parameter of 77 = 10~^Ao. 



scattering also leads to the formation of additional sub- 
gap states bound by multiple Andreev reflection within 
the "pair potential" provided by the suppressed order pa- 
rameter, A||(z). These states appear only near grazing 
incidence and are weakly bound with energies just below 
the continuum edge. 

Sub-gap states do not appear in "^He-A at a specular wall 
since there is no change in phase of the order parameter 
for specular reflection when || z. Thus, all quasiparticle 
states belong to the continuum. This situation changes 
dramatically for a rough surface. Now there are scat- 
tering processes connecting an incident trajectory with a 
reflected trajectory that is at a skew angle, (/> 7^ tt, in the 
xy-plane (see Fig. 0)). For ^He-A the order parameter 
for a trajectory, p = cos^z + sin0(cos0x 4- smcfyy), is 
Aa(p) = z A|| sinSe"^. The change in phase of the or- 
der parameter upon skew scattering leads to strong pair- 
breaking and the formation of sub-gap states. For the 
diffuse scattering the coupling of skew trajectories with 
all possible azimuthal angles generates a band of states 
which fill the sub-gap spectrum as shown in Fig. for 
several incident trajectories. 

The self-consistent spectrum calcuated numerically and 
shown in Fig. [7| is well described by the spectrum ob- 
tained by calculating the retarded propagator for a con- 
stant order parameter, A|[, everywhere in the film. The 
spectrum is then determined entirely by the changes 
in the order parameter induced by diffuse scattering. 
The transport equation can be solved analytically with 
Ovchinnikov's boundary condition for the Riccati ampli- 
tudes. For the diagonal part of the quasiclassical propa- 
gator we obtain. 
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FIG. 7: The density of states in the A-phase for a film of 
thickness D = 4^0 at temperature T — O.OlTc. The excita- 
tion energy is scaled in units of Ao(T), the bulk value of A|| 
at temperature T. 
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where = A(0)^ +£m with A{9) = A|| sin 6*. We calcu- the real energy axis, ie,„ e + iO+. The result for the 
late the retarded propagator by analytic continuation to local density of states is then, 



' \-^{l + t.n^[2uD/vfCos9]) V cos [2^i^M cos 0] 

+e(A(^)^ - e^) , ( ^"^M2^.-I^)Mcosg] _ \ 

1 _ (1 _ tanh^ [2LoD/vf cos 9]) V cosh [2ujD/vf cos 9] J 



where tu is now 

= Vk'-A(0)2|. (31) 

Equation H3U|) shows both the band of the sub-gap states 
generated by diffuse scattering and the Tomasch oscilla- 
tions for energies above the continuum. The positions 
of the maxima are determined by the condition for con- 
structive interference of particle- and hole- like excitations 
with energies above the gap, | A|| sin 6*1 , reflecting from the 
specular surface, 

2./e2- A2 sin2 6' 

D = mT , n = 0,l,2.... (32) 

V f cos 9 



The spectral weight of the Tomasch oscillations also de- 
pends on distance from the surface. Some peaks are sup- 
pressed at special positions in the film due to spatial oscil- 
lations of the particle- hole interference amplitudes. This 
suppression of spectral weight is most visible for angles 
close to 9 = tt/2 near the the free surface. For example, 
the density of states for 9 = OAtt at z = 19/20D shows 
that every second peak is suppressed. 

The sub-gap states are bound to the surface on the 
length scale set by coherence length, ^o, and decay expo- 
nentially into the bulk. However, the situation is differ- 
ent for thin films; the bound states extend over the entire 
width of the film. Figure |S1 shows the total DOS averaged 
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over the film, N(£) = / dz/D J df^p/47rN(p, 2; e). The 
spectrum is gapless over the entire energy range, e < Aq, 
and is finite at e = 0. The inset to Fig. |S1 shows the 
evolution of the total DOS as a function of temperature. 
The gapless states fill the spectrum e < Aq completely 
as T ^ r«"". 

At low temperatures, T — ^ 0, the DOS is insensitive 
to temperature, and the value of the DOS at e = 0, 
< N(0) < N/, persists above e = 0. If we decrease the 
film thickness the sub-gap states fill the gap and the tran- 
sition to the normal state will occur when this process is 
complete. As shown in Fig. |S1 the density of states is 
almost equal to that for the normal state over the whole 
energy range for D = O.S^o- The A to normal transition 
occurs for a slightly smaller film thickness. 

As a consequence of the gapless spectrum the thermo- 
dynamic properties of the films of '^He-A will be very 
different from those of bulk '^He. For example, the low- 
temperature behavior of the specific heat vanishes ex- 
ponentially in bulk ^He-B for T < Ab, C^""'(T) 
{TjTf/^ exp(-Ao/r), while for bulk ^He-A the nodal 
excitations have zero energy at isolated points on the 
Fermi surface. The bulk density of states vanishes at the 
Fermi level as N^""'(£ ^ 0) ~ e^, and the specific heat 
exhibits a power law, C^"'" - {T/T^)^. 

The specific heat of films of superfiuid •^He-A is expected 
to have a different power law behavior at low tempera- 
tures. The density of states is finite and nearly constant 
in the low-energy range above the Fermi level. As a result 
the specific heat will have the linear temperature depen- 
dence as T — > 0, just as for normal '^He, except that the 
Sommerfeld coefficient is reduced in the superfiuid film 
by the ratio, N(0)/N/. 
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THERMODYNAMIC PROPERTIES 
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FIG. 8; The total DOS averaged over film for several film 
thicknesses for T <C T^''". The inset shows the gapless exci- 
tations filling the gap for t = T/T"'" 1. 



To compute the thermodynamic properties of superfiuid 
■^Hc films we need a free energy functional formulated in 
terms of the quasiclassical propagator and self-energies. 
Such a functional has been derived starting from the gen- 
eral Luttinger-Ward functional, formulated in terms of 
the full Green's function and self-energy, by eliminating 
the high-energy, short- wavelength intermediate states, 
and thus computing only corrections to the ground-state 
energy to leading order in the small expansion param- 
eters of Fermi liquid theory. The conceptual formula- 
tion of this problem is discussed in detail by Rainer and 
Serenei^ The formulation of a quasiclassical free-energy 
functional for inhomogcneous equilibrium states is simi- 
lar, but their are additional technical approaches to incor- 
porating inhomogeneities of the order parameter iSSi^2i2i 
Our approach is similar to that of Ref . l26| and is outlined 
in App. m 

We start from the quasiclassical free energy in the weak- 
coupling limit expressed in terms of the quasiclassical 
propagator, 5(p,R;£m) and order parameter, A(p,R), 
derived in Eq. (|ljl6() . 



Af}[5, A] 



rSp'(A5) 



dASp'(A5A). (33) 



The symbol, Sp', denotes the sum over relevant variables: 
the volume of the '^He film, position on the Fermi surface, 
Matsubara energies and a trace over spin and particle- 
hole degrees of freedom. 



Sp'(. 



^^ETr4(...). (34) 



There is an additional integration over the variable cou- 
pling parameter A in Eq. (j^Hfl involving an auxiliary 
propagator, which is the solution of the quasiclassical 
transport equation in Eq. Isj, but with the self-energy 
scaled by A: A ^ A a = AA. The transport equation 
for g\ is not solved self-consistently, but with a single 
integration for each value of A. Thus, g\ is a function 
of the exact order parameter in the film. This procedure 
and the application of boundary conditions for comput- 
ing the auxiliary propagator and the quasiclassical free 
energy functional are also explained in App. ^ 
Equation (|33|l . when evaluated with the self-consistently 
determined propagator and order parameter in the film, 
gives the difference of the thermodynamic potential, 
Af2 = f2s — r^N, from which the change in entropy and 
specific heat can be calculated, 



AS'(r) 



' dT 



aC{T) 



-T- 



(35) 



The normal-state free energy for ■^He of volume V is 
given by n^{T) = Ej, - V{^j„T^), where is the 

ground-state energy for normal "^He and 7^ = ^f-N^ 
is the normal-state Sommerfeld coefficient. 
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FIG. 9: Superfluid correction to the thermodynamic poten- 
tial vs. reduced temperature for several films of superfluid 
^He-A on an rough substrate. The inset shows the reduction 
in the entropy of the superfluid film. 

The reduction of the free energy below the normal-state 
value represents the gain in energy due to the formation 
of a condensate of pairs in the film. The free energy 
of •^He-A films in the limit of diffuse scattering by the 
substrate is shown in Fig. |51for several film thicknesses. 
The reduction in the free energy is given by the Ginzburg- 
Landau form, Ail (x —(1 — T/T"'")^, for temperatures 
just below the superfluid transition temperature, T^''", 
of the film. At low temperatures the gapless excitations 
dominate the thermodynamics. The density of states at 
the Fermi energy is non-zero and approximately constant 
at low energies. As a result the low-temperature limit 
for the free-energy of the superfluid state decreases as 
ils—Es = —V (i7sT^), where 7s < 7n is the Sommerfcld 
coefficient for the low-energy excitations of the superfluid 
film, and Es is the T = condensation energy. 

From the numerical results shown in Fig. ^wc can calcu- 
late temperature dependence of the entropy and specific 
heat of the "^He-A film. The results for the entropy are 
shown in the inset of Fig. |51 The linear temperature 
dependence of the entropy resulting from the gapless ex- 
citations is clearly visible. 

Numerical calculations of the specific heat are shown in 
Fig. ^1 These results show the decrease of the heat 
capacity jump at T^''" with decreasing film thickness, as 
well as the linear temperature dependence of the specific 
heat resulting from the gapless excitations (see inset of 
Fig. O. This behavior for ■^He-A films is in sharp con- 
trast to the low temperature heat capacity of bulk '^He-A, 
which varies as Cs ^ T^. 

Results for the heat capacity jump, AC(r"'")/7ArT"''", 
and the Sommerfcld coefficient, 7s, are summarized in 
Fig. ^2 as a function of the film thickness, D. The Som- 
merfcld coefficient, 7s, was calculated by two indepen- 
dent methods. We calculated 7s directly by numerically 
differentiating the temperature dependence of the super- 
fluid free energy. We can also relate the Sommerfcld co- 




FIG. 10: The specific heat of a superfluid ^He-A film as a 
function of reduced temperature for several film thicknesses. 
The inset shows the ratio, Cs{T)/C,i{T). 



efficient directly to the density of states at the Fermi 
energy, N(0). Thus, 7s/7n = N(0)/N/. The first calcula- 
tion is carried out entirely in the Matsubara formalism, 
while the calculation of the DOS at the Fermi level is 
obtained by the solving for the retarded quasiclassical 
propagator on the real energy axis. Both results agree 
and are shown in Fig. ^2 ^^nd give us confidence in our 
numerical calculations for the propagators, free energy, 
entropy and heat capacity. 




FIG. 11: Specific heat jump, AC(rr'")/7NTc"'" at Tf" 
and the ratio of the low-temperature Sommerfeld coefficient, 
7s /7n, for the superfluid fllm as a function of fllm thickness, 
D. For comparison, AC{Tc)/ji,Tc = 1.19 for bulk ^He-A in 
the weak-coupling limit. 



VI. CONCLUSION 

We have calculated the thermodynamic properties of 
thin films of superfluid '^Hc in the weak-coupling limit. 
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expected to be applicable to films at zero pressure. We 
calculated the phase diagram for the supcrfluid film, in- 
cluding the AB transition, the suppression of the super- 
fluid transition temperature, suppression of the order pa- 
rameter, the quasiparticle density of states and thermo- 
dynamic potential. Our analysis, based on the quasi- 
classical method, shows a spectrum for superfluid films 
with gapless excitations formed by the combination of 
reflection by a rough substrate and Andreev scattering 
induced by changes in the order parameter along classi- 
cal trajectories of quasiparticles. The gapless excitation 
spectrum depends on the film thickness and dominates 
the low temperature thermodynamic potential, entropy 
and specific heat. 
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APPENDIX A: DIFFUSE SCATTERING 

The boundary condition for the quasiclassical propaga- 
tor at a rough surface is obtained by solving the transport 
equation, 

[iemT3 - Si„,p(p, R; £,„) - A(p, Il),g{p, R; £„)] 

+ivf ■ Vg{p,R;e,n) = 0, (Al) 

in the dirty layer and matching this solution to the qua- 
siclassical propagator in the superfiuid. In the limit 
of strong disorder within the impurity layer. |Si„,p| 3> 
\ern\, |a|. Thus, deep in the impurity layer, 

[e..„p(p,R;e™),g(p,R;e™)] =0, (A2) 

with the impurity self-energy evaluated in the Born ap- 
proximation, 

e..„p f ^r-i(p,p')5(p',R;^™) , (A3) 

where T^^(p, p') is the rate for quasiparticles to scatter 
from p ^ p' on the Fermi surface. These equations 
are solved by an isotropic propagator, (7tdl (em), that is 
normalized to g^DL ~ — tt^I. This propagator is not the 
normal-state propagator for the isolated normal metal 
because the proximity coupling to the superfiuid layer 
produces a 'rotation' of gN(ern) 5TDL(£m) in particle- 
hole space. To fix this 'rotation' we include the leading 
corrections to Eq. (|A2|) due to spatial variations of the 
propagator in both the TDL and the supcrfluid film and 
match the solutions at the interface. In the TDL the 
transport equation is 

ivf ■ V3(p,R;£,„) = [Si„p(p,R;e„0,g(p,R;e„i)] . (A4) 



Equations l|A3IA4|l are solved by expanding the prop- 
agator in a basis of Nambu matrices. For superfiuid 
•^He films in zero field the basis is limited to matrices 
in 2 X 2 particle-hole space, with the spin degrees of free- 
dom fixed. Thus, three linearly independent matrices, 
{^1,52,33}, arc required (the identity matrix drops out 
of Eq. (|A4(I . These matrices satisfy the algebraic rela- 
tions of the Pauli matrices, 

= -27r2 6,j , [9t,9j]^ = -2TTe^jkgk ■ (A5) 

We choose 33 ~ ^TDL and express the propagator in the 
TDL as, 

g(p,R;em) = B+(p,R)g+(e„) + S_(p, R) g_(e„J 

+B3(p,R)?3(£m)(A6) 

where 'g± = (pi ± 192) /V2. The linear differential 
equations for {i33(p, R), i3+(p, R), i3_(p, R)} are eas- 
ily solved with Ovchinnikov's model of forward scat- 
tering, T^^(p,p') = 4t~^Pzp'^ for pzp'^ > 0, otherwise 
'''"^(PiP') — 0. Thus, quasiparticles enter the dirty 
layer, scatter forward towards the specular wall, and af- 
ter reflection diffuse out of the TDL. The limit: d ^ 0, 
VfT —^ 0, Vfr/d — > corresponds to diffuse scattering by 
the impurity layer. 

The propagator in the TDL is matched to the propaga- 
tor, 'g{p,0; Em); in the supcrfluid at the interface to the 
TDL. We use the same basis to express 

g(p, 0;£,„) = 5TDL(£,n) + C+{p,0)g+{em) 

+ C_(p,0)5-(£m). (A7) 

The coefficients of this expansion satisfy the following 
relations obtained by Ovchinnikovf24 

C+(p,0) = 0, p, < 0, and (A8) 

/ ^|p,|C+(p,0) = 0, 

C_ (p, 0) = , > , and 

/ ^|p,|C_(p,0) = 0. 

The propagator deep in the dirty layer is also related to 
the physical propagator at the boundary, 

|Pz|?(p,0,em) . (A9) 

Boundary conditions, IA7IA9I can be written in a com- 
pact form using the commutation relations ljA5p . 

?(P, 0;£m) - gTDL(£m) = 

„ . gTDL(£m) , g(p,0;£m) ■ 

This condition is solved self-consistently with Eq. (jA9|) 
for Ptdl and g{p, 0; £„). 
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The boundary condition for ^(p, 0; e„i) can be cast into 
a more compact form using the Ricatti representation for 
the propagator, 5. For an outgoing trajectory, > 0, 
Eq. (|A10|I is solved by 

a(p, 0) = "(iTT - 5tdl)"Vtdl ■ (All) 

Thus, integration along an outgoing trajectory should 
start with the value of a given by the value deep in the 
thin dirty layer. The second Ricatti amplitude, a, is 
known at the TDL substrate, since we integrate a along 
a trajectory with > in backward direction. Thus, 
equations IIA9|l and IjAlip . together with a(p,0;em) and 
the Ricatti parametrization Hl()lll|l are iterated until 
they converge to a value for a(p, 0). 
We use the fourth-order Runge-Kutta method to nu- 
merically integrate the Ricatti equations along a clas- 
sical trajectory for the Ricatti amplitudes, a(p, z; Em) 
and a(p, 2;;em). Azimuthal symmetry for scattering in 
the plane of the film allows us to consider trajecto- 
ries defined by 9 and = (Fig. I12|l . The Integra- 
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FIG. 12: Integration along classical trajectories. 

tion procedure is slightly different for the A- and B- 
phases. In case of the A-phasc we know the Riccati am- 
plitudes at the interface with the substrate (Eq. 125(1 : 
a{d, 0; £,„) = a[9_, 0; Em) = 0. For any trajectory we start 
at point 1 and integrate forward along trajectory 1 — S"— 2, 
with specular reflection at S, to obtain the amplitude d. 
To calculate a we integrate along the same trajectory in 
the reverse direction starting at point 2. 
For the B-phase we do not know the initial values of 
the Ricatti amplitudes anywhere. In this case we start 
with an initial guess for the amplitude a, e.g. at the 
point S. Using Eq. 1161 and inversion in the azimuthal 
plane we find a starting value for the amplitude a. Wc 
then integrate from 5 to 2 to obtain a(^, 0; £m), and from 
S" to 2 to obtain 0(6*, 0; Em). We implement the diffuse 
boundary conditions at the point 2 to obtain d(0, 0; Em), 
and then integrate from 2 to S". This gives us an updated 
initial value for the amplitude a (and by symmetry for a). 
The integration procedure is repeated until convergence 
is reached. 



APPENDIX B: FREE ENERGY FUNCTIONAL 

Wc start with the Luttinger-Ward functional for the full 
Nambu Green's function, G, and self-energy, S, 

n[G,% = -isp{£G + ln(-Goi +S)} + $[G], (Bl) 
where 

Sp{...} = r5: jd'R I ^Tr4{...}, (B2) 

I](p,R;e„i) is the self-energy and Gq (p, R; £,„) = 
i^mTs — Co(p)l is the inverse Green's function for non- 
interacting reference system of bare '^He. The stationar- 
ity condition with respect to the Green's function, 

il^=0^£ = £..[G]=2^, (B3) 

relates the functional $[G] to the self-energy via the 
skeleton expansion for the self-energy, while the station- 
arity condition with respect to the self-energy, 

^ = 0-G-i=Go-i-E, (B4) 

gives Dyson's equation for the Nambu propagator. These 
equations provide a starting point for deriving the Fermi- 
liquid theory for superfluid "^He. In particular, the lead- 
ing order expansion of Dyson's equation in the small pa- 
rameters of Fermi liquid theory can be transformed into 
Eilenberger's transport equation (Eq. [Sj for the quasi- 
classical propagator. 

In order to derive a free energy functional of the quasi- 
classical propagator, g, and quasiclassical self-energy, 

e(p, R; Em) = a (s(p/P, R; em) ^ £n) % , (B5) 

we remove the normal-state stationary point (Gn,Sn) 
from the Luttinger-Ward functional (Eq. IB1|) by defin- 
ing AS = S — Sn, AG = G — Gn and introducing the 
subtracted functional, 

An[G, AS] = n[G, S] - 17[Gf,, En] = 

Sp {as G + ln(-G-i + AS) - In(-G-i)} 

+A$[G] , (B6) 

which has as inputs the normal-state propagator, G^ = 
(G(7^ — Sn)^^, and self-energy, E^, rather than the bare 
propagator. The subtracted ^-functional, 

A$[G] = $[G] - $[G«] - isp {s« (G - Gn)} , (B7) 

is confined to the low-energy region of phase space since 
pairing corrections to the normal-state propagator con- 
tribute only in the low-energy region, ksTc <C Ef. The 
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diagrammatic perturbation expansion for A$ can be re- 
organized as an asymptotic expansion in the small pa- 
rameters of Fermi liquid theorji^ that is formally an ex- 
pansion in the number of low-energy propagator lines. 
To convert Eq. (|B6|) to a functional of the quasiclassi- 
cal propagator and self-energy we integrate out the mo- 
mentum dependence normal to the Fermi surface over 
a region of momentum space near the Fermi surface, 
ICpl < ^c- The low-energy self-energy is slowly varying 
fimction of and can be evaluated with p ~ pjp. Thus, 
the term, AE G, in Eq. ljB6p is <^p-integrated to give 

Sp{aSg} Sp'{65} = (B8) 

d^R / ^Tr4 (s(p,R;£„)5(p,R;£,„)) . 

rn 

To integrate the In-functional wc introduce an auxiliary 
functional defined by introducing a variable coupling con- 
stant for the self-energy and <f>-functional: AE — > AE^ = 
AaE and A$ A$a = Aa$. Thus, the auxihary func- 
tional is 

- isp {aEa G + ln(-G-i + AEa) - In(-G-i)} 
+ Aa$[G]. (B9) 

The stationarity conditions with respect to G and AE> 
give a new equation for an auxiliary propagator, G^^ = 
G~^ — AE>. The auxiliarly functional can be ^p- 
integrated after first differenting with respect to the cou- 
pling parameter, then carrying out the ^p-integration to 
obtain, 

^ = -isp' {eg} + isp' {egx} + A$[g] , (BIO) 
where 

1 f-^"" 

5A(p,R;em) = - / dCpT3 GA(p,R;em) , (Bll) 

is the quasiclassical auxiliary propagator. Wc can inte- 
grate Eq. IjBlOp with respect to the coupling constant. 



Since AJIa^o = and Ar2A=i = Ail we obtain the de- 
sired free energy functional in terms of the quasiclassical 
propagator, g, and self-energy, ©, 

1 

An\g, &] = - dX Sp' {6 (5A - ?)} + A<i>[5] . (B12) 
^ Jo 

The stationarity conditions for the subtracted free- 
energy functional reduce to the quasiclassical transport 
equation and self-energy expansion obtained from the 
asymptotic expansion of the <f> functional. 



g^2^^iM. (B13) 

Sg" 



[i£mT3 -e,g]+ ivf ■Vg = 0, 

These equations are supplemented by boundary condi- 
tions for the propagator g which describe the effects of 
scattering by a surface or interface. 
The auxiliary propagator, g\, is a functional of the exact 
quasiclassical self-energy, and is obtained by solving the 
quasiclassical transport equation with & — > As, 



[iSniTs - Xe,gx]+ivf Vg\ = . 



(B14) 



This auxiliary transport equation is solved once (not 
self-consistently) for each value of A with © as a pre- 
determined input function. The diffuse boundary condi- 
tion for 'gx is given by Eq. (jA10|) with 5tdl fixed by the 
self-consitently determined solution of the quasiclassical 
equations and boundary condition for A = 1. 

A further simplification for AQ is possible in the weak- 
coupling limit when the self-energy is purely off-diagonal 
and given by the order parameter, 6 = A. The self- 
consistency equation. 



A = 2 



can be used to evaluate A$[g] — ^ 
suiting free energy reduces to 



3Sp'{Ag}. 



An 



dASp' A [gx--g 



(B15) 
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